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ABSTRACT 

A substantial fraction of the total stellar mass in rich clusters of galaxies resides in a 
diffuse intergalactic component usually referred to as the Intra-Cluster Light (ICL). 
Theoretical models indicate that these intergalactic stars originate mostly from the 
tidal interaction of the cluster galaxies during the assembly history of the cluster, and 
that a significant fraction of these stars could have formed in-situ from the late infall 
of cold metal-poor gas clouds onto the cluster. However, these models also over-predict 
the fraction of stellar mass in the ICL by a substantial margin, something that is still 
not well understood. The models also make predictions about the age distribution 
of the ICL stars, which may provide additional observational constraints. Here we 
present population synthesis models for the ICL of an intermediate redshift (z=0.29) 
X-ray cluster that we have extensively studied in previous papers. The advantage of 
observing intermediate redshift clusters rather than nearby ones is that the former fit 
the field of view of multi-object spectrographs in 8m telescopes and therefore permit 
to encompass most of the ICL with only a few well placed slits. 

In this paper we show that by stacking spectra at different locations within the 
ICL it is possible to reach sufficiently high S/N ratios to fit population synthesis 
models and derive meaningful results. The models provide ages and metallicities for 
the dominant populations at several different locations within the ICL and the BCG 
halo, as well as measures of the kinematics of the stars as a function of distance from 
the BCG. We thus find that the ICL in our cluster is dominated by old metal rich stars, 
at odds with what has been found in nearby clusters where the stars that dominate 
the ICL are old and metal poor. While we see a weak evidence of a young, metal poor, 
component, if real, these young stars would amount to less than 1% of the total ICL 
mass, much less than the up to 30% predicted by the models. 

We propose that the very metal rich (i.e. 2.5 x solar) stars in the ICL of our 
cluster, which comprise 40% of the total mass, originate mostly from the central 
dumb-bell galaxy, while the remaining solar and metal poor stars come from spiral, 
post-starburst (E+A), and metal poor dwarf galaxies. About 16% of the ICL stars are 
old and metal poor. 

Key words: Clusters of galaxies; diffuse intergalactic medium; stellar populations; 
intra-cluster light. 



discover bridges of light between a triplet of adjacent galax- 
ies in tiie Coma cluster that he attributed to stars "escap- 
ing from the triple system entirely" and making their way 
into intergalactic space. Thus, Zwicky not only pioneered 
the field of difi'use intergalactic light in clusters of galaxies, 
but also correctly guessed the origin of that light. 

Still, the field remained dormant for almost 20 years un- 
til it was revived in the 1970's when it became understood 
that the missing mass in rich clusters of galaxies (which inci- 



1 MOTIVATION 

In the 1930's Fritz Zwicky had a series of profound astro- 
nomical insights. One of them was that intergalactic space 
could not be just empty but had to be populated by stars, 
dust, and gas. This motivated him to look for intergalac- 
tic stars and supernovae in a number of rich clusters using 
first the 18- inch and later the 48-inch Schmidt telescopes on 
Palomar mountain ( Zwicky|1952 1. These searches led him to 
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dentally was also first pointed out by Zwicky in 1933), had 
to be diffuse and not concentrated in the cluster galaxies 



themselves ( de Vaucouleurs & de Vaucouleurs ( 1970 1; Mel 
nick, White, & Hoessel|1l9"77| ) ; |Thuan fc Kormendy| ( |1977| )) 



These early photoelectric and photographic investigations 
indicated that the true intergalactic light did not exceed 
~ 25% of the total stellar light of the Coma cluster and, 
therefore, that the material that binds the cluster had a 
very high mass-to-light ratio. It was also understood in the 
70's that the strong X-ray emission that emanated from rich 
galaxy clusters was produced by thermal Brehmstrahlung 
and, therefore, that the gravitational potential wells of rich 
clusters were filled with large amounts of hot gas. But even 
the mass of hot gas was not enough to solve the missing 
mass problem, and the diffuse stellar intergalactic medium 
was again forgotten for the next 20 years. 

A third revival has taken place in the last decade 
through a combination of realistic numerical simulations and 
powerful new observational techniques. Numerical simula- 
tions have shown that, as Zwicky predicated 50 years ago, 
tidal interactions during the merging history of clusters 're- 
lease' large numbers of stars that make their way into the 
intergalactic space carrying with them the memory of the 
evolutionary history of clusters. In principle, therefore, the 
diffuse stellar component - which we now call the intracluster 
light or ICL - contains a fossil record of the star formation 
history of clusters of galaxies. In practice, however, the com- 
parison between observations and cosmological simulations 
yields surprisingly weak constrains. This is partly due to 
the fact that in both simulations and observations it is ex- 
ceedingly difficult to disentangle bona-fide free-floating stars 
from the halos of the brightest cluster galaxies (BCG), and 
partly to the fact that at a given mass there is a manifold 
of possible merging histories (see e.g. Puchwein et al. (20101 
and references therein). Comparison of theory and obser- 
vations, therefore, require representative samples of clusters 
and thus poses formidable observational (and non trivial nu- 
merical) challenges. 

The current state-of-the-art results indicate, for exam- 
ple, that while AGN feedback reduces the discrepancy be- 
tween models and observations, the ICL fractions observed 
in a reasonably representative sample of clusters are still sig- 
nificantly lower than the fractions predicted by the models 



(Puchwein et al. (20101 and references therein). A surprising 



prediction from these models (with or without AGN feed- 
back) is that as much as 30% of the intergalactic stars are 
actually born in 'cold' intergalactic gas clouds stripped off 
galaxies by tidal interactions. Intergalactic star formation 
has been verified observationally in the intergalactic medium 



(see e.g. Sun et al. (20101 and Boquien et al. (20101 for re- 



cent references), but what is unexpected is that so much 
of the ICL would come from these events. In fact, the dis- 
crepancy between ICL fractions in models and observations 
would basically disappear if these stars did not exist. There- 
fore, studies of the stellar populations of the ICL promise to 
provide useful new constrains to the theoretical models. Un- 
til recently, however, it has only been possible to measure 
the colors of the ICL, which is a notoriously difficult ob- 
servational task. These colors appear to be systematically 
bluer than those of the BCG halos, consistent with a sce- 
nario where much of the ICL originates in smaller galaxies 
(Rudick et al.,||2010|). Also, the ACS on HST resolved the 



stellar populations of the ICL in a field of the Virgo cluster 
and found that, as expected, the intergalactic stars span a 
wide range of ages and metallicities, but the dominant pop- 
ulation accounting for 70-80% of the stars is old and metal 



poor (Williams et al. 20071. A similar result was obtained 



by Coccato et al., (2011 1 in the cluster Hydra-I (A1060) us- 
ing long slit spectroscopy Lick indices. In Hydra-I, basically 
100% of the ICL stars appear to be old and metal poor. 

In this paper we present a spectroscopic analysis of 
the stellar populations in the ICL of an X-ray cluster at 
z=0.29, RXJ0054.0-2823, that we extensively studied pho- 
tometrically and spectroscopically ( Toledo et al.|[20TT |. We 
find that the majority of the ICL stars in this cluster are 
old, as is the case for the Virgo cluster, but the metallici- 
ties are significantly higher indicating a different origin that 
may be related to the fact that our cluster has three inter- 
acting giant elliptical galaxies in the center that conform 
a terrific galaxy grinding machine. The paper is organized 
as follows: In Section [2] we recall the data we used in the 
investigation from previous papers; Section [3] presents pop- 
ulation synthesis models of our stacked spectra of the ICL; 
and Section |4] presents the comparison of our observations 
with previous observations and state-of-the-art cosmological 
numerical simulations, the results and presents the conclu- 
sions of this work. Finally, in Section |5] we summarize our 
results and outline future investigations. 



2 THE DATA 

A detailed description of the data-sets used in the present 
investigation has been presented in three previous papers: 



Giraud et al. (2010 20111, and Toledo et al. (20111; hence 



forth Paper I, to which we refer for comprehensive discus- 
sions of the observing and data analysis methodologies. It is 
useful to reproduce here the relevant figures from Paper I to 
set the data in the present context, in particular because, as 
described above, observationally it is extremely difficult, if 
not impossible, to separate the ICL from the extended halos 
of the brightest cluster galaxies. 

Following previous investigations (e.g. Gonzalez, 



Zabludoff, fc Zaritsky] ( |2005[ ); [Gonzalez, Zaritsky, fc Zablud- 



off ( 2007 1 ) we were constrained to use the ICL and the BCG 
together and to rely on the surface brightness profiles to de- 
lineate the most likely boundaries between the two compo- 
nents. These profiles from Paper I are reproduced in Fig- 
ure [l] The photometric error bars, which include residual 
alignment errors between the V and I images, indicate that 
we detect the ICL out to at least 400 kpc from the cluster 
center. Both the surface brightness and the color profiles ex- 
hibit clear features at about r ~ 50 kpc that we identified in 
Paper I as the radius beyond which the diffuse component 
is dominated by free-ffoating ICL stars. Thus, for r < 50kpc 
we see a mix of BGC and ICL stars, while at r ^ 50kpc we 
observe pure ICL light. We have therefore stacked together 
our spectra in these regions to reach S/N ratios suitable for 
studying their stellar populations, which we discuss in the 
following Section. 
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Figure 1. Left Surface brightness profiles in V and I. Right (V-I) color of the BCG+ICL component as a function of radial distance. 
Special care was taken to match the fits of the two bands to avoid spurious differences due to misalignment and rebinning. Thus, the 
photometric error-bars provide accurate estimates of the reality of the various features visible in this plot. Both plots are on an r^/** 
scale. 



3 POPULATION SYNTESIS 

We fitted simple stellar population (ssp) models to our data 
using the code STARLIGHT described in detail by |Cid Fer-| 
nandez et al. (j2005j). We used a base of 39 ages ranging 
from 6 Myr to 11 Gyr and 3 different metallicities Z=0.004, 
Z=0.02, and Z=0.05, (i.e. metal-poor: 0.2 times solar; so- 
lar; and metal-rich: 2.5 times solar) based on the Galaxev 
code ( Bruzual fc Charlot|2003" I, the Padova 1994 evolution- 



ary tracks (Bertelli et al. 19941, and the MILES library 



i jSanchez-Bl^zquez et al.||2006[ ). STARLIGHT calculates the 
best fitting linear combination of ssp model spectra from 
this BC94c library described in detail by |Gonzlez Delgado| 
|fc Cid Fernandes] (2010) in their comprehensive comparison 
between different stellar libraries on the basis of synthesizing 
the integrated spectra of star clusters mostly in the Magel- 
lanic Clouds. We used the COM extinction law ( |Cardelli| 
et al., 111989] ) and a | Chabrier| p"003| ) IMF, although these 
choices do not affect our results in any significant way. 

As a check, we ran a number of models using different 
stellar libraries (or bases) kindly provided by Roberto Cid- 
Fernandez and partly discussed by |Gonzlez Delgado fc Cid| 
Fernandes (20101 in order to test the sensitivity of our re- 



sults to the choice bases. The best fitting stellar populations 
do depend on the choice of stellar bases, but unfortunately 
these bases do not cover the same wavelength ranges, the 
same age ranges, nor the same metallicity ranges. Therefore, 



as discussed by Gonzlez Delgado & Cid Fernandes ( 2010| , it 
is difficult to choose the library best suited for a particular 
set of observations. 

The four libraries we tested (CB07cQ BC94c; VOOc; 
and VOOs, where s stands for Salpeter and c for Chabrier 
IMFs) yield very similar fits to the data (as indicated by 



^ The CB07 base is similar to the BCOO base of |Gonzlez Delgado| 
|fc Cid Fernandes| poTot but using the MILES library instead of 
STELIB. 



the parameter), but quantitatively different stellar pop- 
ulations. However, all libraries yield the same kinematics, 
and as mentioned above, the synthesis results are insensi- 
tive to the choice of IMF. Thus, our choice of libraries was 
mostly dictated by the fact that BC94 (s or c) span a wider 
range of ages and metallicities than the VOO models. (In 
hindsight, the fact that the BC libraries also span a much 
broader range in wavelength makes them more suitable for 
redshifted objects.) Therefore, while the diffuse light in our 
cluster displays a rich mix of stellar populations with any 
of these libraries, the populations differ in detail for the dif- 
ferent libraries, so one must be careful in interpreting the 
results. 

In order to illustrate the way in which we will present 
the spectral synthesis results for the ICL, we applied 
STARLIGHT to our high S/N spectrum of galaxy #210 (cf. 
Figure [sj which was included in the long slit that we used 
to observe the ICL along the northern side of the major 
axis of the brightest cluster galaxy, which turns out to be a 
dumb-bell system (Paper I). Figure [2] shows in the left pan- 
els the observed (black) and synthetic (red) spectra, nor- 
malized to the fluxes at 400-410 nm, and the flt residuals 
calculated subtracting the model fit from the observed (nor- 
malized) spectra. The magenta colors in the residuals panel 
indicate wavelengths that have been masked from the fit due 
to contamination from poor subtraction of atmospheric ab- 
sorption bands; The red colors correspond to data excluded 
by STARLIGHT using & k — a clipping algorithm. Unfortu- 
nately, the green 01 auroral line at 557.7 nm falls on top of 
the metallicity-sensitive G-band at 430.4 nm, which explains 
why this strong feature has typically high residuals in our 
spectra. 

The right panels of the plot show the population analy- 
sis using different colors for the three metallicities. Present- 
day mass fractions on top, and luminosity fractions on the 
bottom. The labels of these figures show the metallicities 
of the base as well as several parameters resulting from the 
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population fits, notably the velocity dispersion (a), and the 
total visual extinction (Av). The average S/N of the spec- 
trum between 450-500 nm is also given together with an 
estimate of the goodness of fit measured by the reduced x^- 

Galajcy #210 is the fifth ranked elliptical galaxy in the 
cluster with a V magnitude of 19.5, which at the distance of 
the cluster, D ~ 1200 Mpc, makes it an L* galaxy. It's stellar 
populations are dominated by old metal rich stars (92%±2% 
with a small but significant contribution of old metal-poor 
stars 7.5 ± 3%, and an insignificant (at the < la level), 
residual at intermediate-ages. The extinction is negligible as 
expected for this high latitude cluster. In fact, the model 
extinction provides an additional criterion for validating the 
results, since a high value would be indicative of errors in 
the flux calibration of the spectra. 

The velocity dispersion uncorrected for instrumental 
broadening is aobs = 385.2 km/s. The instrumental broad- 
ening measured by fitting a Gaussian profile to the 557.7 nm 
night-sky feature is 15. 5A FWHM, which corresponds to a 
velocity dispersion ainst ~ 354 km/s. Thus, the intrinsic ve- 
locity dispersion of Gal#210 is 152 km/s, typical of elliptical 
galaxies of similar luminosity (e.g. Cheng & Krauss ( 2001 1). 



3.1 Modeling errors 

For a variety of reasons obtaining reliable estimates of the 
statistical uncertainties in the stellar populations derived 
from ssp model fits is a notoriously difficult task (see e.g. 



Gonzlez Delgado & Cid Fernandes (20101 for a recent dis 



cussion) . For spectra with accurate fiux calibrations and rea- 
sonable S/N ratios (say S/N> 15), the random errors are 
small compared to the systematic errors associated with the 
choice of stellar libraries, the wavelength range of the fit, 
the masking of certain features, etc. Thus, understandably 
albeit unfortunately, STARLIGHT does not provide formal er- 
rors in the fitted parameters. 

In order to estimate the intrinsic errors in our popula- 
tions, therefore, we ran each model 10 times to compute 
the dispersion due to numerical errors, and we repeated 
the operation masking some wavelength regions where the 
spectra are particularly noisy, or the residuals particularly 
large. Naturally, the errors depend on the quality of the 
spectra ranging from < 2% for Gal210 to ~ 10% for the 
S-shaped Arc (ARCsum) in population ages; but are signif- 
icantly more stable for the pop. metallicies: < 1% in Gal210 
and < 5% in the Arc. The systematics are sensitive to age 
and metallicity and in general are substantially larger, in 
particular for the ages. In the Arc, for example, the variance 
of different masking experiments can reach up to ~ 20% for 
the ages, but is only ~ 7% for the metallicities. 

A fortiori, therefore, the errors have an objective ran- 
dom numerical fitting error, and a somewhat subjective sys- 
tematic component that is estimated by observing the re- 
sponse of the fits to small changes in the input parameters. 
In the paper we will quote the average values with formal 
one-cr errors, and as relevant, we will comment upon the ef- 
fects changes in the input parameters, which we will refer to 
as systematic variances. 



3.2 ICL population synthesis 

Figure [3] reproduced from Paper I illustrates the way the 
long slits were placed on the cluster. Full details of the ob- 
servations are given in that paper. In order to investigate the 
properties of the ICL, we integrated the pixels at the loca- 
tions indicated in Figure |3] The red spots depict the 4 areas 
that we averaged together to obtain the 'pure' ICL spec- 
trum. Thus, the average of these 4 spectra is what we will 
call ICL-spectrum throughout the rest of the paper. We also 
fitted the spectrum of the halo of the BCG (shown by the 
northern green spot in the inset of Figure [sj , and a mixed 
spectrum of the ICL and the BCG halo indicated by the 
southern green spot in the figure. We will call these spectra 
'inner halo' and 'outer halo' respectively. Finally, we used 
two slitlets to measure the two components of the S-shaped 
arc and we call the sum of these two spectra arcSUM. 



3.3 The BCG Halo and the S-shaped arc 

Observationally it is very difficult to disentangle the bona- 
fide intergalactic component (ICL) from the halo of the 
brightest cluster (dumb-bell) galaxy (BCG; not numbered 
but clearly visible in Figure [3| , which if the cluster has a 
significant ICL, can be very extended. Taking advantage of 
the fact that we used long slits, we are able to sample the 
diffuse light at various distances from the BCG and look for 
systematic trends. 

The left panel of Figure [4] shows the spectrum and the 
ssp fit for the inner BCG halo in exactly the same format 
we used for Gal#210 above. The position of the inner halo 
is indicated in the inset of Figure [3] by the green spot just 
above the Arc. Somewhat surprisingly, while close to 85% 
of the stars in the inner halo are old, only about 20 ± 12% 
are metal-rich while 42 ± 1% are old and metal-poor. The 
reminder have solar metallicity and either old or intermedi- 
ate ages. The old metal-poor population remains stable for 
different experiments, but the systematic variances of the 
other populations are rather large. The strong old metal- 
poor component indicates that a substantial fraction of the 
BCG halo is made of stars accreted from ancient low-mass 
galaxies. 

The velocity dispersion of the BCG Halo (150 km/s) 
is typical of the central velocity dispersions of luminous el- 
liptical galaxies, but seems too large for the outer halos of 



these galaxies (e.g. CaroUo et al., (19951 ). However, in cD 
galaxies the velocity dispersions rise towards the edges of the 
halo ( Kormendy &: Djorgowski,|1989 1 . Unfortunately we do 
not have a complete long-slit spectrum of the BCG to check 
whether we observe a steep velocity dispersion profile, such 
as seen in M87 for example, or a rather flat proflle like those 



observed in normal ellipticals ( Doherty et al. ( 2009 1 and ref- 
erences therein). 

The velocity dispersions of all components are shown in 
Table [T] and the radial velocities of the individual compo- 
nents before averaging is given in Paper I. 

The right-hand panel of Figure |4] shows the spectrum 
of the S-shaped arc (arcSUM), which appears to be a spiral 
galaxy in the process of being disrupted by the tidal field 
of the 'grinding machine'. Spectroscopically the galaxy is 
a typical E-fA or 'post-starburst' galaxy. It shows a very 
complex set of populations, although it should be apparent 
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Figure 2. Left panels. The observed spectrum of galaxy #210 (the fifth ranked elliptical in the cluster) compared to the results of 
the population synthesis STARLIGHT fit in red. The lower panel shows the residuals of the fit. The regions colored in magenta correspond 
to telluric features that have been masked from the fit. The residuals in red have been automatically excluded by STARLIGHT using a 
K — cr clipping algorithm. Right panels. The right panels show the results of the stellar population fit for three different metallicities as 
indicated in the Figure. The upper panel shows the fractional contribution of each component by present-day mass, and the lower panel 
by luminosity. Galaxy #210 is a typical elliptical where most of the stars arc old and metal-rich (red), with a small contribution of old 
and metal-poor stars (blue). 



from the figure that the models are particularly sensitive 
to the masking of difTerent features in the spectrum. The 
most stable population is intermediate age (< 7.5 x 10^ yr) 
and metal rich, with 42 ± 10% of the stars. This component 
could come from a galaxy-wide starburst induced by the 
tidal interaction with the cluster core. About half of the 
stars (with significant systematic variance) are metal-poor 
and have old and intermediate ages, possibly coming from 
the outer arms of the disrupted spiral. The spectrum also 
shows a weak metal poor component at ~ IGyr, which may 
or not be real. If real, it may be a component of the ICL 
rather than of the galaxy itself since the sky subtraction for 
the arc spectra was done far away from the ICL and the 
BCG halo (cf. Figure |3|. 

We remark that our best models for the Arc are not 
able to fit the CN feature at ~ 3800A. We have not found 
a satisfactory explanation for this failure; We looked at the 



populations synthesis fits of thousands of SDSS galaxies ( Cid 



Fernandez 2011) and we did not find a single one with such a 



strong CN feature. We also investigated the possibility that 
in part the feature comes from a higher redshift background 
galaxy. If the CN feature were contaminated by the Mgll h 
and k resonance doublet at 280 nm, for example, in order 



to appear at z=0.29, the redshift of the background source 
would have to be 2 ~ 0.75, and the separation between the 
components of the doublet about 9.5A. Indeed the fit resid- 
uals show a clear doublet structure, but with a separation 
of ~ 40A, so it is not Mgll at z=0.75. An alternative that 
fits this doublet could be Call H & K from a foreground 
galaxy, but we do not observe the expected lines at other 
wavelengths. In the end we opted for masking-out the whole 
feature in order to avoid introducing systematic errors in the 
models. 

The velocity dispersion of the Arc (272 km/s; Table 1) 
is high for isolated spiral galaxies, which should not be sur- 
prising since we are observing the galaxy being stretched 
and distorted by the tidal field of the grinding machine, and 
since we must be picking up a substantial fraction of the 
spiral rotation as well. 



3.4 Moving away from the BCG 

We have two more positions along the major axis (as defined 
by the long axis of the dumb-bell BCG): the outer halo po- 
sition indicated by the green spot in the insert of Figure [3] 
and the crossing of the EW and NS slits also shown in that 
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Figure 3. Approximate positions where the various components of the ICL were measured. For the 'pure' ICL we averaged together the 
spectra in the locations indicated by the red spots. The insert on the top-left of the figure shows the image of the cluster with a stretch 
optimized to show the ICL. For the inner halo of the BCG we used the spectrum at the position of the green spot located just above the 
S-shaped arc in the color insert (bottom-right). The outer green spot indicates the mix spectrum between ICL and halo. We also fitted 
the spectrum at the intersection of the NS and EW slits (clearly visible in the figure marked with a B), and the spectrum of the sum 
the two components of the S-shaped arc (arcSUM). The positions of the slitlets used for sky subtraction are also shown in this figure. 
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Figure 4. Left panel: Spectrum and population synthesis fits of the BCG halo at the position indicated by the green spot above the 
Arc in Figure [s] The strong emission feature at ~ 3800A is not real and was masked from the fit. Right panel: Spectrum and fit to 
the sum of the two arc components (N&S) that were covered by two different slitlets are shown in Figure [s] The strong CN feature at 
380-385 nm could not be fitted by any models and therefore was masked-out to avoid introducing systematic errors. 
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figure indicated with the letter B. The left panel of Fig- 
ure [S] shows the ssp fit for the outer halo. The populations 
are dominated by old metal-rich stars (~ 65% of the mass) 
with a significant fraction ~ 10% of old metal-poor stars. 
The solar-abundance components show the largest system- 
atic fluctuations, but the weak population of the same age as 
the strong burst seen in the S-shaped arc (~ 2.5 Gyr) seems 
to be real. At these positions we are seeing a mix between 
ICL and BCG halo, so the solar abundance component is 
most likely an ICL feature resulting from the disruption of 
spiral galaxies that underwent strong starbursts as they ap- 
proached the core of the cluster, while the metal poor com- 
ponent may arise from the outer parts of these disrupted 
spiral galaxies, or from disrupted dwarfs. 

The velocity dispersion at this position (456 km/s; Ta- 
ble 1) is surprisingly high, higher in fact than that of the 
pure ICL and of the cluster as a whole (see below) . We have 
done several tests by fitting only parts of the spectrum, but 
the velocity dispersion remains high. We do not have a good 
explanation for this high velocity, which we believe to be 
real; we can only speculate that it may be due to the in- 
teraction of the third giant elliptical in the center (#493 in 
Figure [3| with the BCG, which could generate high velocity 
tidal streams. It could also be due to a usual distortion of 
the halo of the central dumb-bell galaxy, which affects the 
velocity field in ways that depend on the actual geometry of 
this subsystem. 

At the location of the slit crossing (r ~ 30 kpc), the stel- 
lar population is completely dominated by old stars of, how- 
ever, difi'erent metallicities. The solar and metal-rich pop- 
ulations suffer the highest systematic variance, while the 
metal-poor component remains stable. In the standard as- 
sumption that the mixed populations are the result of the 
assembly history of the cluster, we must conclude that the 
ICL is formed of stars stripped from galaxies of all sizes. 
However, since most stars have solar and over solar metal- 
licities, the 'donors' of ICL were mostly relatively massive 
galaxies. There are traces (in the luminosity plot) of young 
stars, but typically these appear when STARLIGHT tries to 
fit low-mass hot post-AGB stars that are not represented in 
the stellar libraries. 

The velocity dispersion of the diffuse component at this 
distance from the center (r ~ 40 kpc) is a, ~ 375 km/s 
indicating that this region is mostly dynamically decoupled 
from the BGC and consists mainly of stars freely moving in 
the potential well of the cluster. 



3.5 The (pure) ICL spectrum 

In order to obtain a reasonably high S/N sample of the bona- 
fide ICL, which as we saw in Section [3 . 2| corresponds to r 
50 kpc from the BCG, we averaged together 4 regions where 
the sky-subtracted spectra had reasonably good S/N ratios. 
These regions are indicated by the red spots in Figure |3] 
The resulting combined spectrum shown in Figure [6] has 
S/N ~ 27 between 450 and 500 nm, which seems rather 
good considering the low surface brightness levels involved. 

As expected from the slit-crossing spectrum, the 'pure' 
ICL is dominated by old stars (91 ± 3% of the mass). Of 
these, 38±5% are metal-rich and a similar fraction have solar 
metallicities (with similar errors). As before, the most stable 
population is old and metal-poor, whereas the solar compo- 



nents show the largest systematic variances. There are traces 
of an intermediate-age component of solar metallicity with 
~ 8 ± 3% of the mass, but again with large systematic vari- 
ance. The fraction of luminosity plot (bottom-right) indicate 
traces of very young stars (< 10^) years. These could again 
be due to hot post-AGB stars, but the residuals plot (bot- 
tom left) shows weak nebular emission lines of [OII]3727, 
H/3, and [OIII]5007. It is rather difficult to assess the real 
strength of H/3 because of the strong underlying absorption 
feature, but the line is clearly detected. 

The presence of H/3 is important to understand the ori- 
gin of the nebular lines; they could arise from the interac- 
tion of material ejected from the galaxies with the hot X- 
ray emitting gas in the cluster, or they could depict in-situ 
star formation in the ICL as predicted by the cosmological 
models. The nebular component could also be due to inter- 
galactic planetary nebulae, which are known to be plentiful 



in other clusters (e.g. Doherty et al. (20091 and references 
therein). However, the [OII]/[OIII] line ratios, albeit very 
noisy, resemble more the spectrum of a metal-rich HII re- 
gion than that of a FN. 



4 DISCUSSION 

Table [l] summarizes the kinematical data from the popu- 
lation fits. We see a steep increase in the velocity disper- 
sion with radius qualitatively similar to what is observed by 



Coccato et al., (20111 in Hydra-I. The ID radial velocity 



dispersion of the ICL is 408 km/s with an estimated error 
of ±20 km/s. The velocity dispersion of the BCG halo is 
consistent with the values observed in other giant ellipticals 
and also in the central regions of massive ellipticals. Thus, at 
that position the light is completely dominated by the BCG, 
which seems to be at odds with the population synthesis 
models that show that about 45% of the halo population is 
composed of old metal poor stars, whereas for an elliptical 
galaxy we expect the population to be dominated by metal 
rich stars (as is the case of Galaxy#210 for example). This 
indicates that the BCG halo has a significant component of 
stars stripped from dwarf galaxies. We have discussed above 
the large velocity dispersion of the outer halo, for which we 
do not have a good explanation. 

By definition the ICL stars should be freely floating 
in the gravitational potential of the cluster, and therefore 
should essentially have a similar velocity dispersion to that 
of the cluster as a whole. Interestingly, our cluster has a 
bimodal radial velocity distribution (Paper I), which yields 
two different values for the velocity dispersion depending on 
whether the strongest peak (that contains 90% of the galax- 
ies) represents the cluster, or the double feature is simply 
a statistical fluke (we have redshifts for 95 galaxies in the 
cluster; Paper I). In the first case we get ad — 328 km/s 
whereas in the second case, that is including all galaxies, 
we find ffci = 496 km/s. Our value for the ICL lies exactly 
between these two extremes, confirming that indeed that at 
r 50pc we are seeing the bona-fide ICL population. 

The stellar populations of the ICL are more of a puzzle; 
The state-of-the-art cosmological simulations of IPuchwein] 
et al. (20101 provide detailed predictions for the star forma- 
tion history of the ICL for a model cluster of 10^^h~^MQ at 
z = 0, which is reasonably close to the mass of our cluster 
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ouler.halo EC03 Zbmm out.v9 slit crass BC03 Zbmm out v9 




Figure 5. Left. Spectrum of the outer halo of the BCG as indicated by the lower green spot in Figure [s] Right. Spectrum of the 
position where the southern N-S slit and the E-W slit cross, marked B in Figure |3] The strong unidentified emission-like feature near 
3800A (probably due to an instrumental fluke or improper sky subtraction) was masked from the synthesis as shown in the residuals. 
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Wavelength log Age (yr) 

Figure 6. Stacked spectrum of 4 different regions of the ICL as marked in Figurejs] The subtraction of the synthesis fit allows to unveil 
faint but clearly visible nebular emission lines of [OII]3727, H/3, and [OIII]5007. Although the residual spectrum is too noisy to justify a 
detailed analysis of the nebular parameters, the line ratios are at least consistent with those of metal rich HII regions ionized by young 
massive stars. There are also hints of the presence of such young massive stars in the population synthesis models. 
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Table 1. Velocity dispersions of the studied components 



Component 


R 


O'obs 






kpc 


km/s 


km/s 


Galaxy #210 




385.2 


152 


BCG Inner Halo 


10 


384.3 


150 


S-shaped Arc 


20 


446.5 


272 


BCG Outer Halo 


30 


577.3 


456 


Slit Crossing 


40 


515.9 


375 


ICL 


> 50 


540.0 


408 



(7 X 1O^^M0; Paper I). In particular, their Figure 9 plots the 
distribution of the number of star particles formed as a func- 
tion of redshift that can be directly compared to the results 
of our population synthesis models. The simulations (with 
AGN feedback) predict that about 25% of the ICL stars 
formed more recently than z ~ 1.5. In other words, close 
to 25% of the ICL stars in our cluster should be younger 
than about 6 Gyr. Our observations show that 8 ±3% of the 
stars were formed at t < 6 Gyr, most of which have solar 
abundance. Again the systematic variance of this fraction 
is rather high, but systematic fluctuations never exceed the 
fraction of 8% for these 'young' stars. Moreover, the cosmo- 
logical simulations predict that close to 30% of all stars in 
the ICL are in fact born in situ, in small metal-poor gas 
clouds distributed throughout the clusters. These simula- 
tions also predict that the vast majority of these clouds fell 
into the cluster very recently. Indeed we seem to see traces 
of a hot, young {t < 10 Myr) stellar population, but, if real, 
this population accounts for significantly less that 1% of the 
total stellar mass in the ICL. 

It is interesting to compare our results with the popula- 
tions of the ICL in the Virgo cluster revealed by HST/ACS 
photometry of resolved stars ( Williams et al.||2007 l and by 



the nature of the Arc, which in Faure et al. (20071 we as- 



spectroscopy in Hydra-I (Coccato et al., 20111. The Virgo 



photometry is somewhat difficult to interpret because it de- 
pends on what isochrones they use, and in the way they treat 



AGB stars. Nevertheless, Williams et al. (20071 state that 



in their best fitting models the stellar populations are domi- 
nated by old metal poor stars, which constitute 70% to 80% 
of the total ICL population. There is also evidence in their 
photometry for a younger metal-rich component comprising 
between 20% and 30% of the mass. 

Using Lick indices to interpret their data, [Coccato et al.,| 



( 2011 1 find that essentially 100% of the ICL stars in Hydra-I 



are old and metal poor. In contrast, the ICL in our cluster is 
dominated by old metal rich stars that together with those of 
solar abundance comprise 75 ± 7% of the ICL mass. Thus, 
just like in Virgo and Hydra-I, the ICL population of our 
cluster is dominated by old stars, but in sharp contrast with 
these clusters, less than 20% of the ICL stars in our cluster 
are metal poor. 

It is tempting to speculate that the intermediate age 
stars of solar abundance are the remnants of E-l-A galaxies 
where the A component is the consequence of a tidally in- 
duced, galaxy-wide starburst. In fact, the S-shaped arc in 
the center of the cluster appears to be the remnant of one 
or two spiral galaxies caught in the act of being triturated 
by the grinding machine (Paper I). Our extensive study of 
the cluster was in fact originally aimed at understanding 



cribed to the gravitational lensing of a pair of background 
galaxies. But Figure [4] shows that the spectrum of the Arc 
actually resembles an E-l-A ('post-starburst') galaxy at the 
redshift of the cluster, and where ~ 50% of the stars formed 
in a galaxy-wide starburst 2.5-4 Gyr ago. We stress, how- 
ever, that we do not fully understand the spectrum of the 
S-shaped arc as discussed above. 



5 CONCLUSIONS AND FUTURE WORK 

By stacking our deep spectra of the ICL of an X-ray clus- 
ter a,t z — 0.29 we have been able to obtain a combined 
spectrum with sufficiently high S/N ratio to decompose the 
stellar populations of the ICL. We find that the resulting 
star formation history for the ICL, as told by the population 
synthesis, is not consistent with the predictions of state-of- 
the-art cosmological simulations. While the overall popula- 
tions are dominated by old stars, consistent with the model 
predictions that the majority of the ICL stars come from 
massive galaxies, there is a significant fraction of ICL stars 
that are old and metal poor, indicating that dwarf galaxies, 
and perhaps the outer regions of spiral galaxies, also con- 
tributed significantly to the ICL. Given the structure of the 
cluster, much of the intergalactic stars probably came from 
the interaction of the three massive elliptical galaxies at the 
center, two of which are in an advanced stage of merging 
revealed by their dumb-bell structure. 

We also see an intermediate age population of solar 
abundance, which could be due to stars stripped off E-l-A 
galaxies. This would be the case if spiral galaxies that ap- 
proach the cluster center undergo galaxy-wide tidally in- 
duced starbursts before being disrupted by the grinding ma- 
chine. The population synthesis models also allow us to re- 
veal faint nebular emission lines in the spectrum of the ICL. 
This suggest the existence of intergalactic HII regions, which 
could be ionized by young massive stars formed in-situ. How- 
ever, the total mass fraction of these stars is <C 1% whereas 
the simulations predict that up to 30% of the stars would 
form in the ICL itself. 

It is also possible that the nebular lines arise from un- 
resolved planetary nebulae, although the [OII]/[OIII] ratio 
does not resemble those of planetaries, but the residual spec- 
trum is too noisy to affirm this with any certainty. Finally, 
the nebular emission could be due to the interaction of cold 
gas ejected from the galaxies with the hot X-ray emitting 



plasma that pervades the cluster (e.g. Ferland et al. ( 2009 1). 

Our observations are in sharp contrast with what is seen 
in nearby clusters such as Virgo or Hydra-I where the vast 
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majority of the ICL stars are old and metal poor. The fact 
that our cluster is at a higher redshift {z ~ 0.29) makes no 
difference since galaxy clusters are not observed to evolve 
significantly between z ^ 0.3 and z = 0. It seems more 
reasonable therefore to ascribe the discrepancy to the pecu- 
liar structure of our cluster for which the vast majority of 
the ICL stars were stripped relatively recently after a ma- 
jor merger that brought the 3 giant elliptical galaxies to the 
cluster center. 

It would be of considerable interest to observe a sample 
of clusters at intermediate redshifts {z ~ 0.15 to z ~ 0.25) to 
compare the composition of their ICL for different configura- 
tions of the close environment of the brightest cluster galaxy. 
Observing intermediate redshift clusters allows us to cover 
a substantial fraction of the ICL within the spectrograph 
slits. Using MOS instruments would have the advantage of 
allowing to measure the velocity dispersion of the cluster at 
the same time as we measure the ICL, as we did for our 
z = 0.29 cluster. 
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